Momentum dependence and nodes of the superconducting gap in iron-pnictides 
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We discuss the structure of the superconducting gap in iron pnictides. In the itinerant electron 
picture, gaps with or without nodes have the extended s— wave (s+) symmetry and emerge within 
the same pairing mechanism, determined by the interplay between intra-pocket repulsion and inter- 
band pair hopping. If the pair hopping is stronger, the system develops an gap without nodes. 
In the opposite case the superconductivity is governed by of the momentum-dependent part of 
the pair-hopping, and an s"*" gap shows nodes on electron Fermi surfaces. We argue that the gap 
without/with nodes emerges in systems with a stronger/ weaker tendency towards a spin order. 
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Introduction. The structure of the superconducting 
(SC) gap is one of the most controversial topics in the 
rapidly growing field of iron-based pnictide superconduc- 
tors. Electronic band configuration of the pnictides as- 
sumes two hole pockets centered at the F point q = (0, 0) 
and two electron pockets centered at (0,7r) and (7r,0) in 
the unfolded Brillouin zone (BZ) to which we will refer in 
this paper. Fig. [T^. Multiple Fermi surfaces (FSs) create 
a number of different possibilities for the gap structuresii 

Majority of theoretical works predict gaps of one sign 
on the hole FSs, and of another sign, on average, on 
the electron FSsj^i^i^i^i^ii^i^ We will refer to such gap 
structure as an extended s"*" state. There is no consensus, 
however, on whether the order parameter has nodes 
on electronic FS. Gaps without nodes have been found in 
the 2-band and 5-band itinerant models^ii^iiiSii^ and in the 
localized spin models;^ On the other hand Graser et al 
found^ an s+ state with nodes in the 5-orbital Hubbard 
model. 

Previous results for the gap structure have been ob- 
tained numerically, and there is a call for a simple ana- 
lytical analysis of the pairing mechanisms and the gap. 
In this paper, we show that diflFerent results for the gap 
structure are not fundamentally conflicting, and that 
both nodal and node-less s+ gaps may emerge in the 
same pairing scenario. 

Electron-phonon interaction is weak in Fe-pnictides^ 
and the pairing likely has an electronic origin. In the 
weak-coupling approach, the most natural candidate for 
a pairing glue is the hopping of electron pairs between 
hole and electron pockets {u-^ term in Fig. [T^) . For small 
pockets, a repulsive pair hopping gives rise to an s+ state 
with the sign change of the gap between hole and elec- 
tron FSs, but without nodes, due to small, in parame- 
ter Ep/W, angular variations of the gaps along the FSs, 
where W is the fermionic bandwidth. This simple sce- 
nario is, however, incomplete as there also exists a repul- 
sive interaction within each pocket {u^ in Fig. [T|). This 
interaction does not cancel out from the equation for the 
gap, and we show the s"^ state without nodes is only 
possible when at low energies (smaller than Ep) the pair- 
hopping term exceeds the intra-pocket repulsion. 

The restriction to energies smaller than Ep is essential 
here because bare intra-pocket repulsion is very likely 
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FIG. 1: (Color online), (a) Hole (center) and electron (edges) 
Fermi surfaces in the unfolded BZ and interactions 114 (intra- 
band repulsion) and U3(q — q') (inter-band hopping). The 
s+ SC gap is nearly uniform on the hole FS but may have 
nodes along the dashed lines near the electron FS. (b) The 
folded BZ. Each corner now hosts two electron FSs with gaps 
Aei,2(<^). (c) Schematic RG flow for 114 and the momentum- 
independent part of U3(q — q') up to energies ~ Ep- The 
couplings at this scale determine the pairing. If renormalized 
Us > U4, (line B), the SC pairing occurs even for us{q — q') = 
const, while for M3 < 114 (line A) the pairing is induced by 
the momentum-dependent part in 113(9 ~ <l')j smd Ae{^p) has 
nodes. 



larger than the pair hopping term, and the latter may 
win only if the two interactions are renormalized by high- 
energy fermions and fiow from their bare values at ener- 
gies of order W to new values at energies ^ Ep, Fig. [1};. 
Because of near- nesting between at least one hole and one 
electron FSs the renormalizations in particle-particle and 
particle-hole channels are equally important in between 
Ep and W, and the analysis of the appropriate parquet 
renormalization group (RG) shows^ that the couplings' 
flow is controlled by a density-density inter-pocket inter- 
action which determines the system's tendency towards 
a spin density wave (SDW) order. When this tendency is 
strong, pair-hopping coupling gets pushed up. The intra- 
pocket repulsion meanwhile decreases under RG as one 
should expect for a repulsive pairing interaction. As a 
result, when the system has a strong tendency towards a 
SDW order, the pair-hopping term becomes larger than 
the intra-pocket repulsion somewhere between W and 
Ep, and the system develops an gap without nodes. 
We discuss here what happens for weaker renormal- 
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ization of the coupUng constants, when the pair-hopping 
interaction remains smaller than the intra-pocket repul- 
sion at energies of order Ep. Our key message is that, 
despite strong repulsion, the system still develops a pair- 
ing instability due to momentum-dependent part of the 
pair-hopping term. This instability leads to the s+ gap 
with nodes on the two electron FS and occurs in the 
clean limit at an arbitrary weak pair-hopping. The rea- 
son why the pairing occurs despite strong repulsion is 
quite generic and is related to the fact that the gap that 
oscillates along the electron FS is insensitive to the intra- 
pocket repulsioniii The gap at the hole FS is still affected 
by the intra-pocket repulsion, but it turns out that the 
elimination of the repulsion from just the electron FS is 
a sufficient condition for the SC instability, still driven 
by the pair-hopping term between electron and hole FSs. 

The gap equation. We label fcrmions near different 
FSs as / (hole) and c (electron). For simplicity, we con- 
sider the case when the gaps Ah{q) are equal on the 
two hole FSs, and the gaps Ae(qr) along the two elec- 
tron FSs transform into each other under rotations by 
7r/2, and reflections in all symmetry planes. The ex- 
tension to the case of non-equal gaps on two hole FSs 
is straightforward. For a model of itinerant electrons, 
we consider the interactions within each pocket 114(9 ^ 
q')(ct clqCq'C_q/+/^/lq/q'/_q/) , and bctwecn the pock- 

ets U3{q — q')c\jC^^qfq'f-qi + h.c (we borrowed the nota- 
tions from Ref. [Si). Functions ^3(9 — q') and Ui{q — q') 
have projections onto one-dimensional Aig (s-wave) and 
Big (d— wave) and two-dimensional E (p— wave) rep- 
resentations of the crystal symmetry. We focus on 
the Aig component uf^^{q,q') oc Ya^^ (q)YAjg (q') with 
'^Aigiq) = 1 + b{cosqx + cos qy) {b is a constant). For 
small hole and electron pockets, YA^^iq) ^ 1 + 2b on 
the hole FSs and YA^^iq) = 1 — {bkp/2) cos2ip on the 
electrons FSs, see Fig. [1^. Then ^4(9 — q') w U4, and 
Mg (q, q') « Mg + ^/2u3 cos 2ip (ip corresponds to q or q', 
whichever is on the electron FS). We also studied a Big 
gap, but found that it emerges at a lower T^. 

We assume that the couplings M3, U3 and W4, are al- 
ready renormalized by fermions with energies between W 
and Ep, see Fig. [1]:, and consider the system behavior 
below Ep^ when the pairing channel is decoupled from 
the SDW channel, and the SC instability problem can 
be treated within the BCS approximation. The set of 
the coupled BCS equations for A/j(q) and Ae(q) is ob- 
tained in the standard manner. On the hole FSs Aft,(q) 
is approximately a constant A/j, while on electron FSs 

Ae(q) = Ae -I- Ae(cOSgx -I- COSqy) « Ae -|- Ae'\/2cOS2iy9 

with Ag — ~k'^Ke/2^/2. The coupled linearized equa- 
tions for three gaps A/i, Ag, and Ag are (we take all FSs 
as cylindrical) 



A,, = -u^LAh - u^LAe 

Ae = -U^LAe - U^LAh] 

Ae = -u^LAh, 



U^LAe 



(la) 
(lb) 
(Ic) 



where L = ln(1.13A/Tc), A ~ E-p is the upper energy 
cutoff. From Eqs.[T]we obtain the equation for the critical 
temperature Tc in the form 



(1 + UiLf - ujl^ = u^(l + u4L)L'^ 



(2) 



It is invariant with respect to sign change of U3 and U3. 

The gap that emerges at Tc has generally nonzero Ah, 
Ae, and Ag. Their ratios are 

Ae "3.-, , ,^ Ae U3L 

= — (I + U4L); -— = — — — (3) 

Ae U3 A;, (1 + ^4^) 

where L is a solution of Eq. ([2]). Note that Ae and Ah 
have opposite signs if itg > 0. When the ratio Ae/Ae < 
I/V2, the s"*" state has no nodes, otherwise there are 
nodes on the electron FS. 

Equation ([2]) is a cubic equation in L and can be ana- 
lyzed for arbitrary S ~ U3/U3 and 7 = U4/U3. For 6 — 0, 
it gives L = L() = l/(u3 — U4) = 1/(^3(1 — 7)), and a 
solution for finite Tc exists only when 7 < 1. The gap 
satisfies Ae = — A^ and Ae = 0. 

Consider now 6 0. A naive expectation would be a 
smooth monotonic evolution of the gap with increasing 
S. One indeed finds a smooth evolution, but only for 
7 < 1, when superconductivity exists even at 6 = 0. For 
such 7, L gradually decreases with 6, and the oscillating 
component Ae continuously increases. The nodes on the 
electron FS appear above some critical Scr{j), see Fig. [2] 
We found Scril < 1) « l/\/2 - 7/V^ and 6cr{l) « 
y/8/9{l - 7) for 1 - 7 < 1. 

For 7 > 1, we found a new behavior: the supercon- 
ductivity develops for any S ^ 0. Indeed, the r.h.s of 
Eq. ^ scales as 1*3(7^ - 1)L^ at L > 1, while the 
l.h.s. scales as u^jS^L^. Comparing the two, we find 
L = (7^ — l)/(u3i5^7), i.e. Tc remains finite (but expo- 
nentially small for (5^0): 



Tc = 1.13Aexp 



- 1 



U37(52 



(4) 



The oscillating component of the gap Ae now dominates 
and exceeds both Ah and Ae by 1/S. The reason why 
Tc is non-zero even when the intra-pocket repulsion is 
stronger than the pair hopping (i.e., when 7 ^ 1) is the 
absence of the 1*4 Ae term in Eq. (jlcp because the angu- 
lar integral of Ac{(p) vanishes along the electron FS. In 
other words, oscillations of the gap along the electron FS 
allow a system to avoid strong intra-pocket repulsion and 
develop SC order. 

We see that there is a qualitative change in the sys- 
tem behavior near 7 w 1 and at S ^ 1. In this region 
Ae/Ae = [^2(5^ + (1 - 7)2 + (7 - 1)]/S is nearly discon- 
tinuous, evolving from ©((5) at 7 < 1 to \/2 at 7 = 1, and 
to 0{1/S) at 7 > 1. The ratio Ae/Ah « -1 at (5 <c 1 
and 7 < 1, and decreases as Ad Ah ~ — 1/7 for 7 1. 

On a more careful look, we find two sub-regimes for 
7 > 1, when the gap has nodes. One is the regime of small 
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FIG. 2; (Color online). The phase diagram in dimensionless 
interactions 7 = 114/143 and 5 = ug/ita. The nodeless 
state exists in the lower left corner below the line 5cr(7), when 
intra-band repulsion U4 and momentum-dependent inter-band 
hopping uz are small. Dashed lines show the contours of fixed 
ratio Ae/Ae = 2, 3, 4, 5. The narrow solid line shows 
the position of a minimum of Ae/Ae in the region where 
the gap has nodes. Below this line, the gap develops nodes 
primarily to reduce the effect of intra-pocket repulsion. 
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FIG. 3: (Color online). Top: dependence of the critical tem- 
perature, Tc, in a clean system, on the parameter 7 = M4/U3 
for three values of (5 = 0.1; 0.2; 0.3. For small 5, the critical 
temperature drops fast as 7 approaches 7 = 1. For larger 5, 
the critical temperature remains relatively large even above 
7 = 1. Bottom: contours of constant M3L = M3lnl.l3A/rc = 
4, 8, 16, 32 in the (7, S) plane. 



(5, where the gap is fully adjusted to minimize the effect 
of repulsive U4. In this regime, Ag ^ A^, A/j, i.e., the 
gap is essentially cos 2lp on electron FS and much smaller 
in magnitude on the hole FS. The other is the regime of 
large 5, where the ratio Ag/Ag is again large, but not be- 
cause the system tends to avoid M4 but rather because the 
momentum-dependent part of the effective pair-hopping 
becomes the dominant interaction. In between, Ag/Ae 
passes through minimum, at (5 = ^min(7)- Near 7=1, 
(5min(7) ~ V7 - 1: whilc for large 7, (5min(7) ~ W^l^)l- 
The ratio lS.e//S.h monotonically decreases with increas- 
ing 5 and is not affected by the change of the physics 
upon crossing of 5min(7)- 

The generic phase diagram in the (7, 5) plane is pre- 
sented in Fig. [21 The bolder sohd line is 5cr{l) above 
which s+ gap acquires nodes, and the narrower solid line 
is '5min(7), where the ratio Ag/Ag is at minimum. 

Role of spin fluctuations. Above we considered 

only a direct four-fermion pairing interaction (terms to 
first order in The antisymmetrized pairing vertex 

created from this interaction has equal charge and spin 
components. Beyond first order, such vertex contain non- 
equal spin and charge components, and the spin compo- 
nent describes the interaction mediated by collective spin 
fluctuationsi^ Such interaction is enhanced near an SDW 
instability and, like M3, is attractive for a plus-minus s'^ 
gap. Accordingly, it extends the boundary of an s+ state 
without nodes towards larger 7 Still, unless the system 
is right at the SDW instability point, a SC with a plus- 
minus gap is not realized for sufficiently large U4, and, 
as before, one needs a momentum-dependent part of the 
pair-hopping (or of a spin-fluctuation exchange) to get 



Tg 7^ 0. This, in turn, gives rise to a gap with nodes on 
the electron FS. 

Role of impurities. In the presence of potential im- 
purities, the scattering of electrons results in mixing of 
electron states with different momenta. This mixing does 
not affect superconductors with isotropic order parame- 
ter. The plus-minus s+ SC state at 7 < 1 and (S ^ 1 is 
suppressed by inter-pocket impurity scattering r.n.)^iiiii^ 
but is not sensitive to a much stronger intra-pocket im- 
purity scattering Fg. We now show that Fq strongly re- 
duces the critical temperature Tc(ro) at 7 > 1, when 
superconductivity develops through an anisotropic com- 
ponent of the order parameter on the electron Fermi sur- 
faces. To evaluate the effect of Fq on Tg, we derived 
the self-consistency equations within the Born approxi- 
mation. We found that the linearized gap equations can 
be presented in the form of Eqs. ([T|), but with Eq. (fTa|) 
for the hole FS gap replaced by 

A,, = -u4LAh - usLAe - w3£(Fo)Ag, 

and ip{x) is the digamma function. The equation for 
Tg(Fo) becomes 

(1 + UiLf - ulL^ = ul{l + UiL)LC{To) . (5) 

Both the l.h.s and the r.h.s of Eq. ^ now scale as 
at Tg — > because C{To) tends to a constant at Tg = 0. 
As a result, the solution of Eq. ([5]) does not exist for 
sufficiently small values of tig, i.e., intra-pocket impurity 
scattering is pair-breaking for superconductivity driven 
by momentum-dependent pair hopping. For 7 > 1 and 
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(5 < 1, Tc(ro) vanishes when Tq = Tc/1.13, where is 
the transition temperature in clean samples. 

Application to the pnictides. The electron and hole 
pockets in the pnictides are rather small in size, and it 
is likely that S cx Ep/W is small Our analysis shows 
that for small 6, two situations are possible. When the 
RG renormalizations of U3.4 arc strong enough, due to a 
tendency to a nesting-driven SDW order, 7 = u^/u^ is 
smaller than 1 at energies of order Ep, and the system 
develops a plus-minus gap without nodes and near-equal 
magnitudes along nested hole and electron FS. In this 
situation, the leading instability of undoped system is 
towards an SDW order, so one should expect that it be- 
comes a superconductor only after spin order is destroyed 
by a finite doping which acts against nestingi^. We be- 
lieve this scenario should work for FeAs 1111 and 122 
materials which do show an SDW order at small doping 
and superconductivity at larger dopings. 

If the tendency towards SDW order is less strong and 
M4 — U3 remains positive down to Ep (i.e., 7 > 1), then 
the pairing is predominantly determined by the angle- 
dependent component of the pair-hopping term. In this 
situation, Tc is smaller, and the pairing gap has nodes on 
the electron FS. This is consistent with what is currently 
known about FeP 1111 material LaOFeP which does 
not display a SDW order and has a small Tc ~ 5K. 

This may explain diverse reports of the measured tem- 
perature dependence of the penetration length AA(T) 
The exponential T dependence of AA(T) in SmFeAsO^ 
is consistent with the nodcless Aig gap; linear in T behav- 
ior in LaOFeP— is consistent with line nodes. Measured 
AA(T) in K or Co— doped BaFe2As2 and in F— doped 
La{Nd)FeAsCM scales roughly as T^. The latter behav- 
ior is difficult to obtain from the momentum-dependence 
of the 5+ gap alone, and likely originated from the effect 
of impurities: either inter-pocket impurity scattering F^r 
in an s+ superconductor without nodesr^ or intra-pockct 
scattering Fq for an .s+ SC with nodesj^ 

Comparison with other theories. Our results are con- 
sistent with the numerical RG studies by Wang et al^ 
who found an extended s— wave state, but with a sub- 
stantial modulation of the gap along the two electron 
FS. These authors found that the modulation increases 
with doping, what is in line with the idea that the up- 
turn renormalization of 713, driven by a tendency towards 



a SDW order, decreases with doping. Our results and 
the reasoning are also consistent with RPA-based studies 
of the 5-orbital Hubbard model by Grasier et ali^ They 
considered the case of near-equal intra-orbital and inter- 
orbital Hubbard interactions, which in our notations cor- 
responds to ii3 Si W4. They found that the gaps along 
electron FS have nodes and Ae{(p) ~ Ae(l -I- rcos2ip) 
with r » 2.3 — 2.5. We also found the gap with nodes at 
U3 = U4, and our results are r = 2 at (5 = 0''", and r ~ 2.3 
at J = 0.27. We recall that a small value of 5 is expected 
for small sizes of the pockets because (5 cx EpjW . We 
re-iterate that the gap with no nodes appears only when 
W3 flows to a higher value than U4 under the parquet RG, 
due to non-ladder renormalizations originating from the 
mixing with the SDW channel. Such renormalizations 
arc not included into the RPA formalism. 

Summary. Wc showed here that the gap structure in 
iron-pnictidcs is sensitive to the interplay between intra- 
pockct repulsion and the pair hopping at energies 0{Ep). 
If the pair hopping is larger, the system develops a gap 
without nodes but with different signs of the gap on 
electron and hole FSs. If the intra-pocket repulsion is 
larger, the pairing is still possible, but is now governed 
by the momentum-dependent part of the pair-hopping. 
Such pairing yields an extended s— wave gap with nodes 
on the electronic FS, which allows the system to avoid 
a strong intra-pocket repulsion. We argue that the ex- 
tended s— wave channel always wins over a d— wave chan- 
nel and that the gap with nodes is more likely for the sys- 
tems with less developed tendency towards a SDW order. 
This nodal state is, however, affected by intra-pocket im- 
purity scattering. 
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Note added. A day before this paper was submitted 
to arXiv, T. Maier et al posted arXiv 0903.5216 with 
their analysis of the gap anisotropy, based on numerical 
analysis of a 5 band model. They named the momentum 
dependence of the interaction and the need to overcome a 
repulsion within hole and electron pockets as the reasons 
for the gap anisotropy and the nodes. These are the same 
reasons that wc found in our analytical study. 
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